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This paper presents a package of nine programs for the 
TI-59 calculator. This package was developed as a solution 
to two problems. One problem involved expanding and modify- 
ing an existing set of programs; and a second problem involved 
developing five distribution approximating programs. The 
solution to these problems represents a package with consider- 
able capability in computing confidence intervals, performing 
hypothesis tests and approximating distribution values. The 
distribution approximations include inverse CDF values for 
the Normal, Chi-square, Student's t and F distributions, 
which allow the computation of confidence intervals without 
using tables. 

The TI-59 proved to be a useful tool in solving these 
problems and demonstrated the capability of hand-held programma- 
ble calculators. The comprehensive set of user guides included 
in this programming package provides even the inexperienced 
user with a step-by-step introduction to this capability. 
Additionally, the methods used in preparing this programming 
package are directly applicable to other calculators or 
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I. INTRODUCTION 



The purpose of this paper is to trace the development of 
a package of nine programs for use in the TI-59 calculator. 

Two of these programs compute confidence intervals for either 
the one -population or two -population situations. Another two 
programs perform hypothesis testing, again, for either one- 
or two-population situations. The remaining five programs 
generate approximate distribution values for the Normal, 
Binomial/Multinomial, Chi-square, Student's t and F 
dis tributions . 

The basis for this programming effort was a set of four 
TI-59 programs written by Professor P.. W. Zehna for his 
personal use and later used in the classroom. Professor 
Zehna ' s programs also computed confidence intervals and 
performed hypothesis testing. However, his programs were 
not completely user-friendly, especially in terms of user 
guides; and they were dependent on obtaining some percentile 
values from standard tables. This paper then presents a 
significant expansion in the scope of these early programs. 

The main thrust of this expansion includes a simplified and 
standardized set of programs and user guides, while eliminating 
the dependence on distribution tables. 

This package of nine programs was designed for two types 
of users: the student, who might be asked to quickly solve 

several very different problems in succession; and the working 
analyst, whose main concern is one specific problem which 
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requires great accuracy. Of particular importance to both 
types of users is a detailed set of user guides (Appendix A) 
which includes sample problems. Many of these sample problems 
have been solved in two ways; the first way involves using 
program-generated percentiles; the other way requires input 
of tabled percentiles. While it is not absolutely necessary 
to look up tabled values when using the confidence interval 
programs, that option has been included, and should be used 
when increased accuracy is desired. Except for some cases 
involving small degrees of freedom, however, the accuracy of 
these first two programs is quite good as can be seen in the 
sample problems of Appendix A. When tables are not available, 
the percentile values generated by the five distributional 
programs is outstanding (Appendix B) and can be used in 
lieu of tabled values. One very convenient feature of the 
distribution programs is the ability to provide values normally 
available only by interpolation in the standard tables. It 
should be noted that the methods used to generate approxima- 
tions in the confidence interval program is slightly less 
accurate than those used in the distribution programs. This 
difference is due to the limited number of program steps 
available in the first two programs. Except for the two 
hypothesis testing programs, there is no requirement to use 
the applied statistics module in the Tl-59. However, this 
module is required when performing hypothesis testing because 
of the large program size and the need to provide a signifi- 
cance level with each test. Thus, it was possible to provide 
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detailed, yet simple, user guides to implement theoretically 



correct and accurate programs . 
using either approximations or 
and the analyst can accurately 



By providing the option of 
tabled values, both the student 
solve a variety of problems. 
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II. NATURE OF THE PROBLEM 



The programming effort presented in this paper was 
generated by two major problems. The foremost problem 
involved expanding and modifying an existing set of TI-59 
programs into a more user-friendly package. A second, 
closely related problem was adding the capability for generating 
accurate distribution approximations. The following discussion 
of the solutions to these two problems is a general overview of 
the particular solutions used. The specific methods and theory 
of solution are left for later sections in this paper. 

A. USER-FRIENDLY PROGRAMMING 

User-friendly programming implies programming with the 
user's knowledge, ability, and hardware familiarity in mind. 

In this light, a user-friendly program is a program which 
has significant capability, yet can be used by those with 
only a modest knowledge of either the calculator or the theory 
involved. Applying this definition as a framework for providing 
a user-friendly programming package resulted in four areas of 
effort. These areas are: standardized data entry, standardized 

solution procedures, maximized use of calculator capabilities, 
and improved user guides . 

1 . Standardized Data Entry 

The nine programs in this package all require some 
form of data entry. The data entry schemes for the confidence 
interval and hypothesis testing programs need careful 
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standardization because of the similarity of the data and the 
size of the data sets involved. There are three possible 
types of data in these first four programs. Data may be 
from a one -population sample, a two -population paired sample, 
or a two-population independent sample. Any of these data 
can be easily entered in the form of summary statistics, 
if available; however, raw data requires some standardization 
between programs. The data entry schemes of the first four 
programs use similar data entry subroutines which take 
advantage of the TI-59's data entry sequence. In each of 
these programs, the data entry subroutine is initiated by 
pressing ( D | . The sequence of data entry then differs 
slightly depending on the type of data being entered. All 
of these entry subroutines use a format which requires a 



R/ S to be pressed after each data point entry. This method 



saves one keystroke for each data point compared to the TI-59's 

2 + , entry method. Also, the [r/ S 



two -keys troke , 



2nd 



key is very close to the numerical keyboard, compared to the 
key, thereby eliminating a source of data entry errors 

sequence . 



2nd 



2nd 



2 + 



that frequently occur when using the 

The data entry schemes used for the five distribution 
programs presented no standardization problem. All of these 
programs., except the Multinomial, require only a few parameter 
entries. In the Multinomial program, which can accept as many 
as 35 pairs of parameters, the data entry problem was more 
difficult. Each of the multinomial data points is stored 



separately until computation begins. The data entry 
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sequence automatically repartitions the calculator to make 
room for this potentially large amount of data. 

Data stored directly in registers has been somewhat 
standardized between programs by using the same registers, 
where possible, for similar data. More information concerning 
the contents of data storage registers can be found in 
Appendix A. 

2 . Standardized Solution Procedures 

A considerable effort was made to standardize the 
steps required in each problem solution sequence. In the two 
confidence interval programs (see Appendix A) problems are 
solved in three steps. In the first step, the calculator is 
repartitioned and the data are entered. The second step 
requires entry of either a percentage, in which case the 
program generates an approximate percentile, or a previously 
obtained percentile value. The last step involves the selec- 
tion and initiation of the proper solution subroutine. A 
similar three-step sequence is used in the two hypothesis 
testing programs where the first step includes repartitioning 
and entry of test parameters. The second step involves data 
entry, and the third step, subroutine initiation. 

The five distribution programs contain only a limited 
amount of standardization due to the different nature of the 
programs. The Normal, Chi-square, Student's t, and F programs 
use the HD. DO. CD labels for the same basic functions. 

The label is used in the Normal, Chi-square, and Student's t 

(not F) to generate approximate density values. The [ B [ label 
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is used to generate CDF approximations in all four of these 



inverse CDF approximations. The Binomial/Multinomial program 
shares none of these standardized label uses. 

3 . Maximized Use of Calculator Capabilities 



for the expansion in capability over that achieved by 
Professor Zehna's original programs. Repartitioning allows 
the addition of the program steps necessary for the inverse 
CDF approximating subroutines in the confidence interval 
programs. The added programming space was also used to com- 
pute one confidence interval estimate, for (J u with JJL known, 
which was not available in the original programs. Additionally, 
repartitioning makes possible an F distribution program and a 
Multinomial distribution program. 



larger programs now require three edges of the magnetic 
program cards, which presents a slight inconvenience in 
added loading and storage requirements. Other problems 
associated with repartitioning, such as inadvertent loss of 
program steps and unwanted or improper partitioning, have 
hopefully been eliminated from this programming package by 
extensive validation with sample problems. 

4 . Improved User Guides 

Preparation of improved user guides was a key element 
in making a user-friendly programming package. The complexity 
of the programs involved and the intended use by students 



programs, where similarly, 




The ability of the TI-59 to repartition was the basis 



However, repartitioning is not without its price. The 
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necessitated a departure from the standard TI-S9 program 
record sheet. As can be seen in Appendix A, the improved 
user guides are organized with the user in mind. The general 
outline for each user guide follows this pattern: 

a. Introduction 

b. General Procedures 

c. Specific Procedures 

d. Additional Capabilities 

e. Labels Used 

f. Storage Register Contents 

g. Sample Problems 

The user guides for the five distrubiton programs have been 
combined to take advantage of the relative simplicity in 
using these programs. The answers to each sample problems 
are in a 10-digit format to provide a positive check on 
calculator output when working each problem. The confidence 
interval sample problems are solved in two ways to demonstrate 
the differences between using approximations and tabled 
values for the required percentiles. 

B. ACCURATE DISTRIBUTION APPROXIMATIONS 

The second major problem addressed in this section 
involves generating distribution approximations for both the 
confidence interval programs and the distribution programs. 

The capability to compute accurate distribution approximations 
provides a new dimension to this programming package by 
eliminating the need for standard distribution tables when 
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solving confidence interval problems. A drawback of this 
new capability is the time required to obtain some values 
from the approximating programs. The time required by the 
confidence interval approximations is not nearly as long as 
that for the distribution programs; but then, the quality of 
the approximations is not as good either (see Section III) . 

Appendix B contains comparisons of tabled values with 
both the confidence interval program approximations (Type I) 
and the distribution program approximations (Type II) . The 
actual probability values used in these comparisons were 
obtained using the distribution programs and can be regarded 
as being very close to the actual probability achieved. The 
missing values in the inverse F comparison are due to the 
inability of the Type I approximation to generate inverse F 
values when either of the degrees of freedom parameters is one. 
Also, while only selected approximations are listed in the 
comparison, the inverse CDF approximations are of nearly 
equal quality over the entire range of the appropriate 
function. The Type I approximations displayed in Appendix B 
are generally not as accurate as Type II approximations in 
terms of actual probability achieved. Only the inverse CDF 
approximations for the Chi-square, t, and F distributions are 
presented in Appendix B since all the other approximations 
which are available in the distribution programs duplicate 
table entries. These other approximations include probability 
values, CDF values, and various other distribution values 
(see Appendix A) . The superior quality of these approximations 
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can be attributed to using the same approximating methods 
used in the TI-59 Applied Statistics Module (see Section III). 
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III. THEORY 



This section on theory is presented as a background for 
the solution methods used in the accompanying programs. As 
such, it is not intended to be a primer in statistics. 

Instead, this section should be considered as an intermediate 
level derivation of the specific statistical methods used in 
programming. For a more basic explanation of this material, 
the references cited within each subject area, or equivalent 
texts, should be consulted. 

There are five subject areas discussed in this section. 
First, the theory used for estimating confidence intervals 
in the first two programs is discussed. Next, the hypothesis 
testing theory necessary for programs three and four is dis- 
cussed. And lastly, the methods used in all nine programs 
to obtain approximations to distribution values are discussed. 

A. THEORY FOR ONE -POPULATION CONFIDENCE INTERVAL ESTIMATION 

The derivation of theoretical interval estimates will be 
done in the same order as these estimates appear in the User 
Guide for Program 1 (Appendix A) . Most of these derivations 
use the pivotal -qualtity method to obtain confidence intervals 
(l,u) [Ref. 1: pp . 379-389]. Other methods used here will 

be discussed in slightly more detail, but will still be brief 
compared to the referenced texts. When forming a C.I. with 
a nonsymmetric distribution the interval will represent an 
equal tails solution, where equal tails implies 
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P [X < 1] ~ P [X > u] = (1- Y)/2. This method does not provide 
the shortest C.I. for a given 'y ; however, this method is 
commonly used for its. ease of computation [Ref. 1: p. 382]. 
Regardless of the method used, the resulting C.I. given in 
this section by (l,u) represents the formula used to calculate 
interval estimates. 

1. C.I. for Normal jX with Known [Ref. 2: pp. 77-80] 



(X-yLi)/ (OV^Tn) as the pivotal -quantity , a 100')''% C.I. is 
constructed thus : 



p 

Assuming that X is distributed N ( fj. , (J ) and using 




Substituting z = z ^i + 'y ^/2 = ~ z (l-'Y)/? ancl s i m P li: ^X in g we have 



P X - z(j/^J~n < jU < X + s (J/^~rT = y , or 




2. C.I. for Normal jU with CJ^ Unknown [Ref. 2: p. 80; 



Ref. 3: p. 277; Ref. 1: p. 381] 



Assuming X is distributed and using 

(£-/*)/ (s/^/n), which is distributed t(n-l), as the pivotal- 
quantity, a 100^1 C.I. is constructed thus: 



p t (l-7)/2 (n_l) ^ * fiSt ( 1 + y)/2 J = 7* 



Substituting t = t (i + 'y*) /2 *- n ’ = "‘ t (1 - y) / 2 ^ n " ^ and sim P li_ 
fying we have : 
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or 




3. C.I. for Bernoulli p [Ref. 4: pp . 376-381] 

The pivotal -quantity method does not work for the 
Bernoulli case and another method will be briefly developed 
here. This method starts with two numbers a and b such that 



P 



a < X' < b 




From these limits we have: 



P fx — aj = (1-7* )/2 = P £x S bj , or equivalently, 






na 



= (l-7)/2 = P 



S x i S nb 



Now Y*i i s distributed Binomial (n,p) which 
be explicitly related to the incomplete Beta 
hence to the F distribution to yield: 



for k < n, can 
function and 



fnX - 

u = 

(n - nX) + (nX + 1 ) p ( 1+ 'y )/ 2 ^ v i ’ v 2 ^ 



nX 

1 = — ; : 

(nX + (n - nX + 1 )F( 1 + y y 2 ^ v 2 +2 ‘’ v l~ 2 ^ 

where v^ = (2nX + 2) , and v 0 = (2n - 2nX) 

1 and u form a 1007 % conservative random interval thus: 
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p [1 < p < u] > 7 . 



The outcome (l,u) is a conservative C.I. in the sense that 
for this discrete distribution the confidence that (l,u) 
contains p is at least 100 ')''%. 

4. C.I. for Normal CT 2 with /I Known [Ref. 3: p. 275] 

Assuming that X is distributed NCjLt , CT 2 ) and using 
_ 7 o 

]T( (X^- , which is distributed X ( n ) > as the pivotal- 
quantity, a 1007% equal tails C.I. is constructed thus: 



X 2 d-7)/2 (n) -= Z(v^) 2 /cr 2 -= X 2 (i+7)/2 (n) 



- 7 



Substituting q x = X 2 ( 1 -y ) / 2 (n) ’ and q 2 = X 2 (l*y ) /2 (n) 

and simplifying we have: 



I(X i -,U) 2 /q 2 - a 2 -= I(x r; u) 2 /J = r 



. or 



(l.u> = ( Y. x i 2 ~ n fJL 2 )^/o . 2 » !/<!]_)• 



5 . C.I. for Normal (J 2 with fj. Unknown [Ref. 1 : p. 382; 

Ref. 3: p. 277] 

2 

Assuming that X is distributed N( IJL.cr ) and using 

2 2 2 
(n-l)s / (J , which is distributed X ( n_ l) > as the pivotal- 

quantity, a 1007^ equal tails C.I. is constructed thus: 



p [x 2 (i- 7 )/ 2 (n - 1) •= (n-iw 2 /cr 2 « X 2 (i + 7 )/ 2 {n - 1) ] - 7. 
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Substituting q x = X (1 -*V ) / 2 ^ n_ ’ and q 2 = X (1 + y ) /2 

and simplifying we have: 



(n-l)s 2 /q 2 CT 2 ■< (n-l)s 2 /q 1 



( 1 * u ) . = (n-l)s 2 (l/q 2 , l/q^ ) • 



= 7, or 



6. C.I. for Exponential \ or /I [Ref. 3: p. 279] 

Assuming that are exponential random 

variables with parameter X and using 2 AnX, which is dis- 
tributed X"( 2n )> as tJie pivotal -quantity , a 100y% equal 
tails C.I. is constructed thus: 



X 2 (i-y)/ 2 (2n) / 2nX mg:2 '^ nX X 2 (i+y)/ 2 (2n) 



=r- 



Simplifying we have: 



X 2 (i-y)/2 (2n) / 2nX -A* X 2 (i + y)/2 (2n)/ 2nx 



-X 



or 



(i,u) = ( X 2 (i_y )/ 2 (2n)/2nx » X 2 (i+y )/ 2 ^ 2n V 2n x) • 

The C.I. for the mean time to failure ( M - i/X ) 15 
constructed by inverting the above interval to yield: 



(l,u-) = (2nx/X 2 ( 1 + y j/ 2 (2n) , 2nx/X 2 ( 1 _y y 2 (2n) ) 

|Ref. 4 s p. 362J . 
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B. THEORY FOR TWO-POPULATION CONFIDENCE INTERVAL ESTIMATION 
The confidence interval estimates for two -population 
situations are discussed here in the same order as they 
appear in the User Guide for Program 2 (Appendix A) . All 
of these estimates use the pivotal -quantity method discussed 
earlier [Ref. 1: pp. 379-389]. As in the one -population 

case above, the C.I. given by (l,u) represents the formula 
used to calculate the interval estimate. 

1. C.I. For Bernoulli P^-Py f° r Large m and n [Ref. 2: 
p. 249] 

Large m and n means np x> mp Y ,n(l-p x ), and m(l-py) all 
greater than five. With this condition met and assuming that 
X and Y are independent and normally distributed, 



which is. distributed approximately N(0,1), is used as the 
pivotal -quantity . A IQO'Y % C.I. is constructed thus: 



(X - Y) - (p x - p Y ) 





(X - Y) - (p x - p Y ) 




Substituting c = ->/x(:L -X)/n + Y(l-Y)/m , and simplifying yields 



P X - Y - cz^ 1 + y <- P x ~P Y ■< X - Y + cz 



d +7 )/2 = 7 > or 
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(l,u) = (x-y - cZ( 1 + yy 2 , x-y +cz ( 1 + y )/ 2 ^ • 



2. C.I. for Normal Mx‘ My for X and Y Paired [Ref. 2: 
p. 123] 

Assuming that X and Y are normally distributed and 
letting D = X - Y, and using (D- *'/n) , which is 

distributed t(n-l), as the pivotal -quantity , a 1007% C.I. is 
constructed thus: 



P . 



t d-y)/2 (n ' 1) ^( D -M D )/s D /Vn) < t d+y )/2 (n_1) = *^ 



Substituting - /Z y , and t = t + / 2 Cn-1) = 

' t (l-7)/2 (n_1) anc ^ simplifying yields: 



D - tS j/yn Mx'My D + ts D / V^] = ^ 



, or 



Zif-W 2 /* 

(l,u) = (d - tSjy/ yn, d+tsQ/^/n) , where s D = [^TT) 



and d. = x i - y ± . 

3. C.I. for Normal My with C^ 2 = (Ty 2 = CJ 2 

[Ref. 2: p. 123] 

Assuming that X and Y are independent and using 
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, which is distributed t(m+n-2), 



(* - * ) - 

(l/m + l/n)S 2 

Jr 

as the pivotal -quantity . Here and elsewhere in this paper, 
n and m represent the number of data points in X and Y, 
respectively. A 1007% C.I. is constructed thus: 



(x-y) - (U x -n y ) 

' { i ^ / 9 (m+n-2) < 5 ^ ‘t/ 1 , *y ^ 

U-7)/2 (l/m+l/n)S 2 (1+ 7>/2 



= 7 ' 



Substituting t = t ^ 2 (m+n- 2) = - 1 ^ _<y ^ / 2 (m+n- 2) and 
simplifying yields: 



X - Y - t V(lA+ l/n)S 2 < /i x - ^ Y «=: X - Y + t l/m + l/n)S 2 =7, 



or 



(l,u) = (x-y - t Vo-/ m + l/n) s 2 , x -y + t l/n)s 2 ), 

V ir 



Bxj-x) 2 + pyj^-y) 



- \2 



where s 



( m + n - 2 ) 



4. 



2 2 

C.I. for Normal [JL y with <J X and (Jy Known 

[Ref. 2: p. 123] 

Assuming X and Y are independent and using 



(X - Y) - ( fJi x - j ± Y ) 

V OxV.-i + CTy^/m 



, which is distributed N(0,1), as 



the pivotal-quantity, a 100 7% C.I. is constructed thus 
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„ (x - y) - ( n Y -LL y ) 

Z (l + 7)/2 «< 7 = = = ^j= =j L < Z (l + 7)/2 



Vcr x 2 /n + (Jy 2 /™ 



= 7 - 



Substituting z = z ^ + < y)/ 2 = " z (l-'y )/ 7 anc * simplifying we 
have : 



X - Y - z V CT x 2 /n + C Y 2 /m < /X x ~ /i Y < X - Y + z V Cr x 2 /n + 0^ 2 /m 






or { 1 , u ) = (x - y - z V CT x 2 /n + O^ 2 /™ > x - y + z V C x 2 / n + (7 Y 2 /m) • 



5. C.I. for Normal (J x 2 / Cy 2 [Ref. 4: p. 464] 

Assuming X and Y are independent and using 
(Sy"/ CTy 2 ) / (S x ^/ (Jy“) , which is distributed F(m-l,n-l), as 
the pivotal -quantity , a 100*X% C.I. is constructed thus: 



s Y 2 /<7 Y 2 



F (l-y>/2 (, "' 1 ' n - l) * - 2, *2 * F a + 7)/2 (m - 1 ’ n - 1) 

b x /u x 



= 7 - 



Substituting F ( 1 .y3/ 2 (m-l,n-l] F ( . i + y )/2 (n-i ,m-l) ’ and 



simplifying we have: 



2/o 2 



S x /Sy 



CT X 



(S 2 /S 2 )F 



^(l + T)/2 U - 1 ’ m - lj (7 2 '“X/-Y ^(i + 7)/ 2 



(m-1 , n-1 ) 



= 7 - 
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or (l,u) = ( 



2 / 

s Y /s. 






(s x 2 /sy 2 )F 



(i+y )/ 2 (m-1 * n_l ) )• 



6. C.I. for Exponential Xx^Xy = E Ref * 4: 

p. 466] 

Assuming X and Y are independent and using X Xx /Y Xy> 
which is distributed F(2n,2m), as the pivotal -quantity , a 
100^% C.I. is constructed thus: 



F (l_y )/2^ 2n,2m ^ ^ x Xx/^Xy ■< F ( 1 + y ) / / 2 ^ 2n,2m ^ 



= r- 



Simplifying we have: 



(Y/^)F (l _y } y 2 (2n,2m) <X x /Xy < (Y/x)F (l + y )/ 2 (2n, 2m) 



=r- 



or 



(l.u) = ( ( Y/X)F {l _yy 2 (2n,2m) , (Y/X)F (1 + y y 2 (2n,2m)) 



C. THEORY FOR ONE-POPULATION HYPOTHESIS TESTS 

Three types of hypothesis tests are performed in Program 3; 
these are upper- tailed , lower-tailed and two-tailed tests. The 
general procedure used for all three tests is the same. 
Basically, a test statistic T, which has an assumed distribu- 
tion, is computed with user-supplied data and then compared 

to a critical region defined by T . This T value is deter- 

0 o o 

mined by the type of test, the assumed distribution, and the 

user-supplied OC value. For two-tailed tests there are 

actually two values of T used: these values will be denoted 

o 
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here by and T ? . In performing the test then, if the test 
statistic T is outside the critical region, then we accept 
our original hypothesis H q ; otherwise, we reject H Q and 
accept the alternative H^. 

The following graphs illustrate the three types of tests 
discussed. These graphs represent a probability density 
function, f (T) , where the shaded area under each curve 
represents values from a corresponding cumulative distribution 
function , F (T) . 



Upper-Tailed Test 
f(T) 




Accept H q — Reject H q 



(Critical Region) 
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Lower-Tailed Test 
f(T) 




Two-Tailed Test for Symmetric Distributions (Normal, t ) 

f(T) 




Two-Tailed Test for Unsvmmetric Distributions (Chi-Sauare. F ) 
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An upper-tailed test is performed whenever the user 
enters a +1, choosing the upper- tail alternate hypothesis, 
during step one of the solution. A lower-tailed test requires 
the user to enter a -1 during the solution, and a two-tailed 
test requires a 0 entry. For two-tailed unsymmetric tests 
the relationship between the test statistic and the median 
determines which tail is actually used for the test. If 
the test statistic is greater than the median then the two- 
tailed test is performed using the upper tail T Q value. For 
test statistic values less than the median the test uses the 
lower tail T Q value. However, there is no general agreement 
that a two-tailed test should be performed this way. 

The test statistics used in Program 3 will now be listed 
in the same order in which they appear in the User Guide. 

This listing will also include a reference where more informa- 
tion concerning a particular subject can be found. Additionally, 
the assumed distribution for each test statistic will appear 
with that test statistic. 

2 

1. Test Statistic for Normal {J. Q with (J Known 



(Ref. 1: p. 431] 



T 




, using N (0 , 1) 



CT 



2. Test Statistic for Normal jj[ Q with Unknown 




[Ref. 1: p. 431] 
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T 



using N(0,1) for n > 30, and 
t(n-l) for n < 30 



(x - /I 0 )VT 

s x 

3. Test Statistic for Bernoulli P q [Ref. 2: p. 101] 

For n > 30 
n(p_ - x) 

T = " ~~ » using N( 0, 1 ) 

V n(l - P 0 )P 0 

For n < 30 the Binomial Distribution in the Statistics 
Module is used to directly calculate the appropriate probability 
for comparison with CX . 

7 

4. Test Statistic for Normal (J Q with JJ. Known 

[Ref. 1: p. 432; Ref. 2: p. 104] 

For n < 65 



T = T' = 



IV + n ^ 2 - 2 M£*i 



using X"Cn) 



For n > 65 



5 . 



m _ T' - n 

Test Statistic for Normal 
[Ref. 1: p. 432; Ref. 2, 
For n < 64 



using N(0,1) 

? 

<V with Unknown 
p. 104] 
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(n-l)s v 2 9 

T = T’ = 5 ^- , using (n- 1) 

0 * 2 

For n > 64 

T = — j.H~ 1 A , using N(0,1) 

V 2 (n-l) 

6 . Test Statistic for Exponential Mo ■ !/Xo [ Ref - 3 
p. 279] 

For n < 32 



T = T ' 



_ 2 nx 

' Mo 



using X"( 2n ) 



For n > 32 



T = 




using N(0,1) 



7. Test Statistic for Poisson [Ref. 2: p. 248] 



For n < 30 

T - n X 0 



using ( 2 nx) 
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For n > 30 



T = 



n * - n X 0 

V^To 



using N(0,1) 



D. THEORY FOR TWO -POPULATION HYPOTHESIS TESTS 

The two -population hypothesis tests are performed in the 
same manner as the one -population tests. For an explanation 
of these tests, refer to C above. The test statistics used 
in Program 4 will now be listed in the same order in which 
they appear in the User Guide. This listing will include 
the assumed distribution and applicable references. 

1. Test Statistic for Bernoulli Py = Py [Ref. 2: p. 249] 



T 




_ ini 



+ nx 
(m + n) 



(x - y) 

)\ Amy + 



“)C 



(m + 





using N(0,1] 



2. Test Statistic for Normal = jJL y for X,Y Paired 



[Ref. 2: p. 121] 

T = H’(/n/s^) , using t(n-l) for n < 31, and 

N (0 , 1) for n > 31 
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3. 



Test Statistic for Normal = JJy for X and Y 

2 2 

Independent with (Jy = (Jy“ [Ref. 1: pp . 434-435 ; 

Ref. 2: p. 116] 



X - Y 



T = 



./ni 5 

o\Jn m 



where s = 



'(n-l)s x 2 + (m-ljsy 2 
n + m - 2 



and using t(n + m-2) for n < 32, and 
N(0 , 1) for n > 32 

4. Test Statistic for Normal Mx a h Y for X and Y 

^ 

Independent with (Jy , (Jy Known [Ref. 2: p. 119] 




, using N(0 , 1) 



5. 



Test Statistic for Normal JJy = JIy 

7 7 

Independent with (Jy f (Jy [Ref. 



for X and Y 
2: p. 119] 




, using N(0,1) for df _> 30, and 
t(df) for df < 30, where 
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/ 



df = largest integer in 




- 1-5 



) 



2 2 

6. Test Statistic for Normal = (Jy for X and Y 

Independent [Ref. 2: p. Ill] 



T - 




using F(n-l,m-l) 



7 . Test Statistic for Normal P = 0 [Ref. 1 : pp . 492- 

493; Ref. 2: p. 200] 




using t(n-2) for n < 28, and 
N (0 , 1) for n > 28 



8. Test Statistic for Exponential = Xy ^ or ^ anc * ^ 
Independent [Ref. 4: p. 310] 

T = — , using F(2n,2m) 

y 



E. THEORY FOR METHODS OF APPROXIMATION 



The methods used to approximate distribution 
this package of programs come from many sources, 
complete discussion of the theory involved would 
the scope of this paper. Therefore, while every 



values for 
and any 
be beyond 
method which 
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has been used will be referenced, only those methods unique 
to this programming effort will be discussed in theoretical 
detail. With this in mind, a combination listing and dis- 
cussion will follow which will trace the approximation methods 
used in this paper. The approximations unique to this effort 
are the inverse cumulative distribution functions (CDF) for 
the Chi-Square, Student’s t and F distributions as well as 
the Multinomial approximation. The inverse Normal CDF 
approximation discussed here was used in Professor Zhena's 
programs. All other approximations used in this paper are 
modifications of methods from the TI-59 Statistics Module. 

For a better discussion of the exact methods used in these 
cases, the TI-59 Applied Statistics Manual should be consulted 
[Ref. 7]. 

There are two types of inverse CDF approximations used in 
this paper. The first two programs use a less accurate 
Type I approximation, while the distribution programs use 
a closely related, but more accurate, Type II approximation. 
What follows is a discussion of the methods used for 
approximating the inverse Normal, Chi-Square, Student's t and 
F distrubitons , and the Multinomial distribution. 

1. Inverse Normal CDF Approximation [Ref. 5: p. 933] 

a. Type I Approximation 

This approximation is used in Programs 1 and 2 
for inverse Normal CDF values and as a subroutine for the Chi- 
Square, t and F approximations in those same programs. The 
Type I approximation uses the following set of equations and 
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